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ABSTRACT 

We use N-body simulations to investigate the evolution of the orientation and magnitude of dark 
matter halo angular momentum within the large scale structure since z=3. We look at the evolution of 
the alignment of halo spins with filaments and with each other, as well as the spin parameter, which is 
a measure of the magnitude of angular momentum. It was found that the angular momentum vectors 
of dark matter haloes at high redshift have a weak tendency to be orthogonal to filaments and high 
mass haloes have a stronger orthogonal alignment than low mass haloes. Since z=l, the spins of low 
mass haloes have become weakly aligned parallel to filaments, whereas high mass haloes kept their 
orthogonal alignment. This recent parallel alignment of low mass haloes casts doubt on tidal torque 
theory as the sole mechanism for the build up of angular momentum. We see evidence for bulk flows 
and the broadening of filaments over time in the alignments of halo spin and velocities. We find a 
significant alignment of the spin of neighboring dark matter haloes only at very small separations, 
r < 0.3Mpc/h, which is driven by substructure. A correlation of the spin parameter with halo mass 
is confirmed at high redshift. 



1. INTRODUCTION 

The large scale structure of the universe observed to- 
day has formed by a long history of gravitational collapse, 
gradual accretion and mergers. Through these processes 
a filamentary, sponge-like structure has emerged. The 
distribution of galaxies and their motions provides clues 
on how they formed, and together with galactic angu- 
lar momentum data, the emergence of the intricate large 
scale structure can begin to be explained. 

Before we can determine what spin tells us about the 
formation of large scale structure, the mechanisms of an- 
gular momentum build-up need to be well understood. 
The initial spin of earl y dark matter proto-haloes can be 
predicted analytically ( W hite] |1984[) . however these pre- 
dictions are largely limited to the regime of linear struc- 
ture formation. To track the angular momentum build 
up through more recent cosmic history, N-body simula- 
tions of cold dark matter must be used. These simu- 
lations give full information on the dark matter haloes 
which can be used to form hypothesis on the build up of 
galaxy angular momentum on cosmological scales. How- 
ever, on cosmological scales it is not yet feasible to sim- 
ulate the gas component to track the angu lar momen- 
tum b uild up of galaxies directly (although Hah n et all 
(2010) simulated 100 disk galaxies in a filament, to find 
an alignment of galaxy spin with filaments). 

Hydr odynamical simulations on individual galaxy 
scale s (Ivan den Bosch et aTl 120031 : ISharma fc Steinmetzl 
120051 : iBett et al.ll2010D have shown that the specific an- 
gular momentum of baryons remains close to that of dark 
matter and that the galaxy angular momentum is gen- 
erally about 20° misaligned with the dark matter halo. 
This means that dark matter halo spin is a fairly good 
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proxy for galaxy spin, so some understanding of the spins 
of galaxies may be gleamed from dark matter-only simu- 
lations. The spin of a dark matter halo depends mainly 
on two things; the initial torques driven by the sur- 
rounding landscape at early times, and the accretion and 
merger history of the halo. 

The initial spin of dark matter haloes is given through 
a mechan i sm kn o wn as "tidal to rque theory" , pi o neere d 
bv lHovia (fl94^ . lF^ebIesl (fl969l) and IZelTJovichl <fl970h . 
This theory proposes that the initial spin of a proto-halo 
early in its formation in the linear regime of structure 
formation depends on its shape and the tidal forces ex- 
erted from the surrounding structure, so the spin is de- 
pendent on the local dark matter landscape. The great- 
est effects of tidal torquing happen at the time of turn- 
around, just before the proto-haloes have collapsed to 
virialized objects. A halo that was torqued in this man- 
ner should retain some memory of the tidal field where 
it formed, and this has been confirmed through N-body 
simulations and galaxy catalogues (eg. ILee fc Pen! 120011 : 
iPorciani et ail 120021: ILee fc Erdogdul 120071 ). The cosmic 
web is the manifestation of the tidal field, filaments in 
particular are regular, symmetric morphologies which on 
large scales exhibit a uniform tidal field. Thus it is ex- 
pected that the orientation of halo spin today should 
retain some correlation with the direction of filaments 
and haloes should be aligned with each other over short 
distances. 

Since the epoch of tidal torquing, halo spins have been 
substanti ally influenced by me rgers and accretion. It was 
shown in IBett & Freiikl (|2011D that it is not uncommon 
for the direction of the spin of a halo to completely flip 
over in its lifetime and this phenomenon is caused by mi- 
nor and major mergers and even close halo flybys. Satel- 
lite accretion has been proposed to be the main con- 
tributor of angular momentum and it has been shown 
that by neglecting tidal torques and considering merg- 
ers alone, the distribution of the magnitude of spin can 
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be reproduced (se e lGardner|[200ll: iVitvitska et"aHl2002t 
IMaller et alJ[200l . 

To figure out how accretion has influenced dark mat- 
ter halo spin and what spin can reveal about the forma- 
tion of large scale structure, several authors have inves- 
tigated an alignment of spin with the cosmic web using 
N-body simulations and galaxy catalogues. In simula- 
tions, it has been found that spins are aligned on shells 
around voids, lying preferentially on th e void surface 
(jBrunino et al.l 120071 : ICuesta et al.1 12008). It has been 
shown that spi ns lie preferentially in the plane of sheets 
in simulati ons dNavarro et a'l . 2004) a nd along the axis of 
filaments dFaltenbacher et al. I2002t lAragon-Calvo et al.l 
I2007bt iHahn et al.ll2007ri IZhang et all 12009ft . In obser- 
vations there has been a te ntative detection o f some weak 
correlation with filaments (jJones et~aT1[2010l ) but no sig- 
nificant detection has been found to date. The evolution 
of halo spi n with respect to fi laments and sheets was ex- 
plored by IHahn et all (|2007aT ) who found no change in 
the orientation of spin over cosmic time. 

Since the spins of haloes are aligned with the large 
scale structure, there should be some degree of coher- 
ence between the direction of spin of two neighbouring 
haloes. It is not clear if this alignment is st rong enough to 
be det e cted even in N - body s imulations. iHeavens et all 
(20001 . iPoreiani et al.l (|2002|). [Faltenbacher et al.l (|2002[ ) 



and " iBailin fc Steinmetd (l2005fi " see no strong alignment 
whereas lHatton fc Nininl (|2001| ) do see a significant align- 
ment. In contrast, several claims have bee n made of 
spira l galaxy spin alignments in ob servations (|Pen et al.1 
[20001: ISlosar et al.l 120091: lLe^l201lh . If these alignments 
can be seen in observations but not in dark matter sim- 
ulations then it is a possible indication that the spins 
of the luminous galaxies are not aligned with their dark 
matter haloes. 

As well as the orientation, the magnitude of the spin 
may reveal secrets of the large scale structure. The spin 
parameter is a dimensionless measure of the amount of 
rotation of a dark matter halo and it has been found 
(jLemson fc Kauffmaimlll999tlCervantes-Sodi et alJl2008| ) 
not to depend on co smol ogy or environment. Both 
iKnebe fc Powerl (j2008l ) and lMunoz-Cuartas et all (|2011[ ) 
find a mass dependence of the spin parameter at high 
redshift but not at low redshift. 

Observations of galaxy spin alignments in the large 
scale structure to date have only been through inferred 
galaxy spin o rientations from observ ed disk galaxy shape. 
For example, iLee &: Erdogdul (120071) us e d the Tul l y cat- 
alogue of nearby spirals (|Nilsonl 119741 : lLaubertsI Il982f) 
to infer spin from the axial ratio (to find an alignment 
with the tidal field) and ISlosar et al.l (|2009l ) used the ap- 
parent sense of spiral rotation in the Galaxy Zoo cat- 
alogue. Direct measurements of galaxy rotation have 
been done with integrated field units (IFU) although 
only one galaxy is targeted at a time and it is not fea- 
sible to conduct a survey of large scale structure with 
direct spin measurements. However, a new multi-object 
IFU instrument has been developed which will enable 
a survey of 10 galaxies in a volume limited sample 
([Bland-Hawthorn et alJfcOllUCroom et al.ll2011h . There 
will soon be a huge influx of galaxy spin data, which 
has never been sampled before in such high volumes. In 
order to get the most out of these data and to direct fu- 
ture surveys, the dark matter halo spin must be better 



understood. 

Our paper is organized as follows. First, the method 
are described in Section [2j Here we describe the set of 
simulations used in Section [2 .11 then we discuss the char- 
acteristic mass scale for halo collapse in Section 12.21 and 
the method used for finding features in the large scale 
structure are described in Section 12.31 Theoretical pre- 
dictions from Tidal Torque theory are discussed in Sec- 
tion [3] and the results of alignment of halo spin with fil- 
aments and the alignment of neighbouring haloes' spins 
are presented. Results of the evolution of the spin param- 
eter in are presented in Section [4] Lastly we summarize 
and discuss our results in Section [5] 

2. METHOD 
2.1. N-Body Simulation 

Since any relic alignments of spin with the large scale 
structure are expected to be weak, a large simula- 
tion volume and high resolution are needed. To this 
end, the publ i cally available Millennium simulation of 
ISpringel et all (2005) was used. This simulation is of a 
cubic volume 500Mpc/h on a side containing 216 3 par- 
ticles using the GADGET- 2 code (|Springelll2005l ). This 
gives a particle mass of 8.6 x 10 8 M Q /h. A ACDM cos- 
mology is chosen and the parameters are fi m = 0.25, 
n b = 0.045, fl A = 0.75, h= 0.73, n = 1 and cr s = 9. 

The halo catalogue was built bv lSpringel et al.l (2005) 
by first using the simple friends-of-friends group (FOF) 
finder (|D avis et al]|1985|) to attempt to select structure 
in the particle distribution and then finding the virial- 
ized subhaloes within the FOF groups using SUBFIND 
(jSpringel et al.ll200l . The SUBFIND algorithm first 
identifies subhalo candidates within each FOF halo us- 
ing dark matter density and then removed particles that 
are not gravitationally bound to the subhalo candidate. 
The most massive subhalo typically contains most of the 
mass of the corresponding FOF object, and so can be re- 
garded as the selfbound background halo itself, with the 
remaining subhaloes as its substructure. The halo cata- 
logue used in this paper includes all virialized haloes, in- 
cluding subhaloes, although spin measurements are only 
made on haloes with more than 500 particles in order to 
minimize random effects from outer halo particles. There 
are 184,891 FOF haloes and 213,799 haloes in total. 

For this analysis, a 300Mpc/h section of the full Mil- 
lennium simulation was used. This smaller section was 
chosen so that the resolution of the density field was high 
enough to be able to find features in the large scale struc- 
ture. This was tested using several lOOMpc/h sample 
cubes. As the resolution of the density field was raised 
from 128 3 to 1024 3 voxels, the alignment between halo 
spin and the resulting filaments became stable above a 
certain threshold. For smoothing lengths 2.0, 3.5 and 5.0 
Mpc/h (Gaussian smoothing is used for finding filaments 
on different scales, see Section |2~3)) . the minimum reso- 
lution for stable features is 0.4Mpc/cell. For a grid of 
1024 3 voxels, the maximum box size is 400 Mpc/h. To 
ensure the resolution was more than sufficient, a box of 
size 300 Mpc/h was chosen. For smoothing on 1.0 Mpc/h 
scales, a finer grid must be used and the maximum cell 
size is 0.2 Mpc so a 200 Mpc/h box was used for this scale. 
At smaller scales than 1 Mpc/h the boxsize required is 
too small so there are not enough haloes for useful re- 
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Fig. 1. — Left: The distribution of dark matter haloes in a volume of the simulation where the large scale structure has been dissected 
into its component features. Haloes in blob regions are colored black, filament haloes are dark gray, sheet haloes are light gray and haloes 
in voids are outlined in black. The size of the dots are proportional to the virial radius of the halo and the volume shown is 100 X 100 X 5 
Mpc/h. Right: The volume of the simulation is uniquely classified into features of the large scale structure using the dark matter density 
field. Here the classification of the volume is shown through the shading: blob regions are black, filaments are dark gray, sheets are light 
gray and voids are white. The features have been found on the scale of 2Mpc/h. 

sheets and voids. This analysis uses the curvature of the 
density field to identify each of these features in N-body 
simulations. 

Firstly, the density field is obtained using the 
Delaunay Tessellation Field Estimator (DTFE) 
method using the dark matter halo di stribu - 
tion (see Ivan de Weygaert fc Schaapl (120070 : 
ISchaap fc van de WeveaertJ (|2000fl : ISchaad (|2007l) ). 
The DTFE method can be summarized in three steps, i) 
from the distribution of points the Delaunay tessellation 
is constructed, which is a volume covering division of 
space into mutually distinct Delaunay tetrahedra. A 
Delaunay tetrahedron is defined by the set of four points 
whose circubscribing sphere does not contain any of 
the other points in the generating set. ii) The local 
density at each point is calculated from the volume of 
the Voronoi cells (the dual of the Delaunay tessellation) 
and the mass of the contained halo, iii) The density 
within each Voronoi cell is interpolated, assuming the 
density field varies linearly. The DTFE method is useful 
when looking for geometrical features in the density field 
because it automatically adapts to variations in density 
and geometry. 

The DTFE was carried out with vacuum boundary 
conditions and a buffer region around the box. This 
buffer region was made to be at least as big as the max- 
imum distance between nearest neighbor haloes so that 
no voronoi cells constructed leaked outside the filled re- 
gion. For larger smoothing scales, the buffer was at least 
as big as 2a. For the 2 and 3.5 Mpc scales the buffer was 
7Mpc and for the 5Mpc scale the buffer was f0.5Mpc. 
The buffer region was also used in the smoothing of the 
density field then discarded. 

Smoothing the density field to some scale s is done by 
convolving with a spherically symmetric Gaussian filter, 



suits. The following results display no cosmic variance 
when a different sample of the same size is chosen. There 
are 4,027,242 haloes in our 300 Mpc/h box and 932,96f 
haloes with more than 500 particles from which a reliable 
spin measurement could be made. The haloes in a 5Mpc 
slice through the simulation volume are shown in Figure 

Snapshots are taken at several points throughout the 
simulation. Here we have used the snapshots at redshift 
0, 0.99, 2.07 and 3.06 (rounded to 0, 1, 2, 3). 

2.2. Characteristic Mass 

In structure formation, there is a characteristic mass 
scale for collapse, M*(z). A spherical top- hat perturba- 
tion collapses when its linear overdensity exceeds a value 
of 8 C = 1.686. The variance of linear density fluctuations 
at a given mass scale M is related to the linear power 
spectrum P(k, z) at redshift z by 

i r°° 

a 2 (M,z) = — dkk 2 P(k,z)W 2 H (k,M), (1) 
27r Jo 

where WTn(k, M) is the Fourier transform of a spher- 
ical top-hat window function of comoving size R = 
(3M / 47T/?) 1 / 3 , and p is the comoving mean mass den- 
sity of the universe. At a given redshift, the typical mass 
scale M*(z) to collapse from a la fluctuation is hence 
given by the implicit solution of 



a(M*,z) = 5 C 



(2) 



The calculated values of characteristic mass at redshift 
0, 1, 2 and 3 are 5.89, 0.273, 0.0132, 4x 10~ 5 , respectively 
in units of 10 12 M Q . 

2.3. Quantifying the large scale structure 

Morphological features in large scale structure may be 
classified into four general categories: blobs, filaments, 



As(x) = J dyp(y)G s (x,y). 



(3) 
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Here p{y) is the Fourier transform of the DTFE density 
and the Gaussian filter at scale s is defined by, 



1 
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The curvature of the density field is given by the Hessian 
matrix of second derivatives at each point, 



d 2 Ps {x) 

dXadxp 



(5) 



The second derivatives can be found while simultane- 
ously smoothing the field by making use of an identity of 
the convolution; J^(/*S , ) = fj*fl , = /*gf- Applying 
this to Equation [Sfgives 



This feature finding algorithm uniquely identifies re- 
gions into blob, filament, sheet or void depending only 
on the scale and quality of features required. 

3. ALIGNMENT OF HALO SPIN WITH THE COSMIC WEB 

Halo particles can be loosely bound, following stochas- 
tic paths, but adding up each particles angular momen- 
tum gives the nett effect of a halo spin. Spin is calculated 
by adding up the angular momentum of each particle (i) 
in the halo, simply defined as the cross product of the dis- 
tance of the particle from the halo's centre of mass (r) 
and the particles velocity (v) with respect to the centre 
of mass; 



£- 

i=0 



Yi x rriiVi 



(8) 



d 2 Ps (x) 
dx a dxa 



[ d 2 
= / dyp{y)- — - — G s (x,y). 

J dx a dx/3 



(6) 



Thus, the Hessian of the smoothed density field is simply 
given by the convolution of the DTFE density and the 
second derivative of the Gaussian (the so-called 'Mexican 
Hat wavelet'.) 



H, 



1 
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In order to get a reliable measurement of halo spin, only 
the haloes with more than 500 particles have been in- 
cluded. The unit spin vectors are shown in the top panel 
of Figure [2] but there is no obvious alignments with each 
other or with the large scale structure (as defined by the 
axis of filaments, shown on the bottom panel). 

From tidal torque theory (TTT), the spin of dark mat- 
ter haloes is expected to be correlated with the local tidal 
dyp(y)[(x a - y a )(xp - yp) - S a/3 a 2 }G s (7) field (T = Ty = didj<f>) and the inertia tensor (I = ly). 

During the linear regime (assuming that T and I are un- 
correlated), the first order result from TTT (|Whitdll984l) 
is, 

Ji oc eijkTjiIi k , (9) 

where is the Levi-Civita symbol. In the principle axis 
frame of the tidal tensor, where Ai are the eigenvalues of 
the tidal field, 



The eigenvalues of the Hessian quantify the curvature of 
density at a particular point, in the direction of the cor- 
responding eigenvector. A positive eigenvalue indicates 
that the shape of the density field is concave up and a 
negative is concave down. The density field may now be 
classified uniquely into blob, filament, sheet or void re- 
gions according to the eigenvalues of this Hessian. The 
eigenvalue sign criteria for each region is as follows, 

Blob All negative 

Filament Two negative, one positive 

Sheet Two positive, one negative 

Void All positive 

It can be useful to cl assify every point in to one of these 
features as was done in lZhang et al.l (|2009j) . and an alter- 
na tive approach is to p i ck out o nly the best features like 
in lAragon-Calvo et all ()2007aD . The decomposition of 
volume into features is shown in Figure [T] on the scale of 
2Mpc/h. The filament and sheet morphologies dominate 
the volume, with blob regions taking up the least vol- 
ume. The relative volume fractions do not change much 
over scale. 

Morphological features are defined using only the 
eigenvalues of the Hessian. The direction of the eigenvec- 
tors are also used to assign a directionality to filaments 
and sheets. The direction of the axis of a filament is 
the direction of the positive eigenvalue, and the normal 
direction of a sheet is the direction of the negative eigen- 
value. The features discussed in this paper have been 
found choosing the smoothing scales of 2.0, 3.5 and 5.0 
Mpc/h. These scales have been chosen to mat ch with the 
visual classification of structure at 2Mph/h (jHahn et al.l 
l2007bf) and to explore the scales above that. The comov- 
ing smoothing scales are kept constant for different red- 
shifts in order not to bias the results with preconceived 
assumptions about filament formation. 
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A3 < A 2 < Ai so A 3 — Ai is the largest coefficient, making 
J2 the largest component of J so that spin is preferen- 
tially aligned with the second eigenvector of the tidal 
field. The cosmic web is a manifestation of the poten- 
tial (f>, related by the Poisson equation, \/ 2 cj) = AirGp (x). 
Our definition of a filament (having two negative eigen- 
vectors of the Hessian of density) translates into a region 
where there are two positive eigenvectors of the tidal ten- 
sor. The second eigenvector of the tidal field points in 
a direction orthogonal to the filament (the minor axis of 
the tidal field is the axis of the filament) and so we expect 
that halo spin should point in a direction orthogonal to 
the axis of the filament. 

The result from TTT in Equation [9] assumes that T 
and I are completely un correlated, which has been shown 
to be not always true (|Lee fc Pen! [2000; Por ciani et al.1 
2002). If there is some correlation, the preferred direction 
of halo spins discussed above may be a small effect. The 
alignment would also be greatly affected by merger and 
accretion events that have happened during nonlinear 
structure growth. 

An expression for the relation between the unit spin 
vector ( J) and the unit traceless tidal field (T) was pro- 
posed in [LeE&Pin| HH HHI): 
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Fig. 2. — The direction of dark matter halo spin vectors (top), 
velocity vectors (middle) and filament axis (bottom). The veloci- 
ties show a coherent flow along filament axis whereas spin vectors 
are much more random and not obviously aligned. Shown is a slice 
of the simulation 100 X 100 X 5 Mpc/h and all vectors have been 
normalized to have the same length. 



where c £ [0, 3/5] is the correlation parameter to measure 
the strength of the intrinsic spin-shear alignment with 
the nonlinear modifications taken into account. When 
c=0 it corresponds to the case when nonlinear effects 
have completely broken down initial spin-shear correla- 
tions and when c=3/5 it is the ideal case when I is inde- 
p endent of T . 

iLee et all (|2005D derived an expression using Equation 
[TOl for the Probability Density Function (PDF) of the 
orientations of the galaxy spin vectors relative to the tidal 
spin tensors: 

1 3 

P(cos a, cos cos 9) = — TT(1 + c - 3cA- )~ 1/2 x 



cos 2 (3 



+ 



cos 6 



1 + c - 3cA? 1 + c - 3cAi 1 + c - 3cA| 



-3/2 



(11) 



Where \ are the eigenvalues of T and a, (3 and 8 are the 
angles between the unit spin vector and the major, in- 
termediate and minor axis of the tidal field, respectively. 

To quantify the preferred alignment of halo spins or- 
thogonal to filament axis, we calculate P(cos 9) which is 
the PDF of the cosine of the angle between spin axis and 
the minor axis of the tidal field which defines the axis 
of filaments. Filament regions are defined as having two 
positive and one negative eigenvector. They also must 
satisfy the traceless condition of J^. Aj = as well as the 

unit magnitude condition of Yli tf — 1- Therefore the 
eigenvalues in filament regions can be approximated by 
Ai = A 2 = 1/V6 and A3 = — 2/V6. Using these values in 
Equation [TT] gives 



c I 1 cos 

2 



-3/2 



P(co S 0) = (l-c) A /l + - 1 

( 12 ) 

If halo spins are oriented completely randomly then c = 
and the PDF is flat. If halo spins are preferentially or- 
thogonal to filaments then c > and the function in- 
creases with cos 6. Although tidal torque theory restricts 
c to positive values, other effects could be in play that 
cause halo spins to be aligned parallel with filaments, 
which would cause a negative value of c. 

3.1. Alignment of halo spin and velocity with filaments 

The alignment between a filament and the spin of the 
haloes that make it up is simply given by the cosine of 
the angle 9 between the two vectors and the absolute 
magnitude is taken because the filament is only defined 
by an axis, not a particular direction. The distribution 
of I cos f?| for all haloes in filaments at redshift and 3 is 
shown in Figure [3] where the number of haloes in each bin 
of I cos 9\ is normalized to make the area under the graph 
unity. The shape of this distribution can be quantified 
in two ways; the median value or by fitting a function to 
the curve. 

Since the distributions shown in Figure [3] are clearly 
non-Gaussian, the median rather than the mean would 
be the more useful s tatist i c (although the mean was us ed 
bveg. iZhang et all (|2009f ); IArag6n-Calvo et all (|2007bD ). 
The standard error of the median was found by bootstrap 
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Fig. 3. — The distribution of alignments of the spin with the 
axis of filaments of all dark matter haloes. The solid lines are for 
z=0 (these are all the haloes in the red region in Figure [5j and 
the dashed lines are for z=3 (the haloes in the green region in 
Figure [5}- For z=0, the overabundance of haloes with high values 
of | cos0| (where 8 is the angle between halo spin and filament axis) 
indicates that haloes are preferentially aligned parallel to filaments. 
For z=3 there is an alignment of spins orthogonal to the axis of 
fila ments. The red lines are fits from theory (the PDF in Equation 
[HI where c is the best fit value, c = -0.035 ± 0.004 for z=0 and 
c = 0.129 ± 0.009 for z=3). The flat dashed line is the expected 
distribution for random halo spin orientations. 



resampling and finding the standard deviation of the re- 
sampled medians. The distributions can also be fitted to 
the probability density function of Equation to find 
the correlation parameter c of the intrinsic spin-shear 
alignment which characterizes the shape of the distribu- 
tion. The fit was done using a Markov chain Monte Carlo 
and two examples of such a fit is shown as the red lines 
in Figure 21 

These above two methods are compared in Figure U 
for some example points (haloes in mass bins at z=0, 
scale= 2.0 Mpc which are the same mass bins as the red 
line in Figure [5]) . There is a one-to-one correlation of the 
two parameters so either could be used. We have chosen 
to use the correlation parameter c in this paper since it 
is theoretically motivated by TTT. 

The value of c indicates the strength of the alignment 
of halo spins with the orientation of filaments, and also 
the intrinsic alignment of spin with the tidal field. If 
the haloes generally have spins parallel to filament axis 
c is negative, conversely, if the halo spin are generally 
orthogonal to filament axis then c will be positive. The 
error of c is the standard deviation of the value which 
maximizes the likelihood of the fit of the PDF to the dis- 
tribution. From the value of c found for all the haloes at 
z=0 (c = -0.035±0.004) and for z=3 (c = 0.129±0.009), 
the general trend is that haloes are aligned orthogonal to 
filaments at high redshift and aligned parallel at low red- 
shift. 

The alignment of halo spin vectors with filaments is 
shown in Figure [5] The alignment distribution has been 
fitted to find c for haloes in bins of mass and for haloes 
at different redshifts. For all smoothing scales, it can be 
seen that at z=0 the alignment is weakly parallel (neg- 
ative c) for low mass halos in filaments (mass less than 
about M* = 5.89 x 10 12 M Q ) and orthogonal (positive 
c) for high mass haloes. This is illustrated in Figure [SJ 
At higher redshifts the alignment becomes more orthog- 
onal for all halo masses. There are less haloes in the 
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Fig. 4. — A comparison of two ways to quantify the degree 
of alignment of halo spin with the large scale structure. The 
data points are for haloes in mass bins at z=0, scale= 2.0 Mpc 
where the same mass bins can be seen in the red line in Figure 151 
(log(M) =11.63 - 12.21, 12.21 - 12.79, 12.79 - 13.37, 13.37 - 13.95, 
13.95 - 14.53). The error bars of c are the la errors of the MCMC 
fit and the error bars of < |cos6| > are the standard error of the 
median. The dashed lines are for random spin orientations. 



high mass bins at high redshift because the high mass 
haloes have not h a d tim e to form yet. The r esult o f 
Faltenbacher et all (120021) ; lAragon-Calvo et all (j2007blk 



Hahn et all (|2007bD and iZhang et all (|2009f ) that halo 



spins generally lie along the axis of filaments is driven 
by the low mass haloes at z=0. This is demonstrated in 
Figure [3] where the alignment distribution for all haloes 
at z=0 is shown. The alignment is preferentially paral- 
lel because of the high number of low mass haloes that 
exhibit parallel alignment. 

The affects of smoothing scale on the halo spin align- 
ment with filaments show something about the forma- 
tion of filaments. For redshift (the red line in Figure 
haloes seem to be best aligned at a large smoothing 
scale while high redshift haloes are best aligned at small 
scales. If an orthogonal alignment is an indicator that a 
halo formed inside a filament topology, then this shows 
that filaments grow in size over time. 

Figure [7] shows the effect of taking into account the 
characteristic mass. Here we can compare haloes be- 
tween redshifts at equivalent stages of collapse. When 
the this is accounted for, almost all the points overlap 
within their errors. This means that haloes at a similar 
stage in their collapse have the same degree of preferen- 
tial alignment with filaments over cosmic time. A halo 
that is just starting to collapse (M = A/*) at redshift 
2 has a similar probability of orthogonal alignment with 
its filament as a halo that is just starting to collapse at 
redshift 1 or 0. However, no assumptions were made 
about the evolving scale of filaments and the smoothing 
scale was kept constant at 2.0Mpc. Even with a con- 
stant scale, this similarity between alignments at differ- 
ent times shows that the build up of spin is closely linked 
with a halo's formation. 

When substructure is discounted by taking the most 
massive subhalo in each FOF group, there is practically 
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Fig. 5. — T he a lignment of dark matter halo spin with filaments over cosmic time. Alignment is characterized by the parameter c of the 
fit of Equation ll2l to the distribution of | cos0|, where positive c indicates orthogonal alignment and negative c indicates parallel alignment. 
The panels show filaments found in different smoothing scales: 1.0 (top left) 2.0 (top right), 3.5 (bottom left) and 5.0 Mpc/h (bottom 
right). At high redshift all spins are orthogonal to filaments but recent times, low mass haloes have a parallel alignment with filaments. 
The dashed line is the expected distribution for random halo spins and the shaded regions are the la errors. The red line is for z=0, yellow 
line is z=l, blue is z=2 and green line is z=3. 



no change in the alignments. 

Although the c parameter was introduced in the con- 
text of spin alignments with the tidal field (manifested by 
filaments in the large scale structure), it can also be used 
as a more general measure of alignment. The distribu- 
tions of | cos^l where 9 is the angle between halo centre 
of mass velocity and filament axis is also well-fit by the 
PDF in Equation [T^l Again, a negative value of c means 
a parallel alignment and a positive value is orthogonal 
alignment. 

All panels of Figure [5] show a parallel alignment which 
is stronger for high mass haloes. This shows streaming of 
haloes of all masses down filaments into massive clusters. 

This streaming can be seen in the velocity vectors of 
haloes in some filaments in the middle panel of Figure 
[21 where vectors are pointed along filaments towards 



clusters. However, some filaments display bulk motions 
where the entire filament is moving towards some attrac- 
tor. To see the extent of these bulk motions, they have 
been subtracted in Figure [9] by subtracting the mass- 
weighted average velocity of haloes by halo mass found 
within the smoothing scale on which the filaments were 
found. When bulk motions are discarded, an orthogonal 
motion remains. The apparent streaming of haloes down 
filaments was wholly caused by bulk motions of entire 
filaments, and this bulk flow is generally along the axis 
of filaments. The relative motions can be seen in Figure 
ITD1 in the alignment of halo velocity with the flow of the 
local bulk motion. (Bulk motions have been subtracted 
from halo velocities here.) Low mass haloes are moving 
slightly orthogonal to the flow and high mass haloes have 
no preferred direction of motion. This reflects how bulk 
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Fig. 6. — The distribution of alignments of halo spin with the 
axis of filaments for low mass (M < M„, solid lines) and high 
mass (M > M*, dashed lines) haloes. These haloes are at z=0 and 
filaments are found using smoothing scale 3.5Mpc/h. This mass 
division corresponds to the red line in the bottom left panel of 
Figure \5\ where the first two points of that figure are the low mass 
haloes and the other points are the high mass haloes. 
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Fig. 7. — The alignment of dark matter halo spin with filaments 
over cosmic time for haloes in bins of halo mass normalized by the 
characteristic mass. The alignment for a 2Mpc scale is shown here. 



motions have been removed: high mass haloes were given 
more weight than low mass haloes and so the residual 
motions of high mass haloes once bulk flow is removed is 
minimal. 

The enlargement of filaments over time that was seen 
in the spin alignments is also visible in the way the bulk 
flows are aligned. The low mass haloes at z=0 (red line 
in Figure [8|) are more strongly aligned at large smoothing 
scales and the low mass haloes at high redshifts are most 
aligned at small smoothing scales. If filaments are chutes 
where haloes are channeled into clusters then these low 
mass haloes are evidence for the growth of the size of fila- 
ments over time. The high mass haloes on the other hand 
are generally less aligned at large smoothing scales for all 
redshifts which is seen as a flattening of the curves. This 
may be due to the inclusion of some cluster haloes when 
the smoothing scale is broadened which would introduce 
random velocities into the sample. 



Although both halo spin and velocity are somewhat 
aligned with filaments, these alignments are not strong 
enough so that there is a significant alignment between 
a halo's spin and velocity. 

3.2. Halo-halo spin alignment 

Tidal torque theory predicts that as well as being 
aligned with the large scale structure, halo spins should 
be aligned with each other. This is usually tested by 
simply taking the average of the dot product of pairs of 
halo spins separated by distance r; 

r,(r) = (\3(x)-3(x + r)\). (13) 

A second quant i ty us ed by iPen et all (|2000D and 
iBailin fc Steinmetzl (|2005D is 



m (r) = (\3(x)-3(x + r)\ 2 



(14) 



These quantities are plotted in the top panels of Figure 
[T2l where at very small halo separations (r < 0.3Mpc/h) 
there seems to be a parallel alignment of halo spins. 

However, both of these quantities rely on taking an 
average over all the halo pairs in each bin of separation. 
The mean is a useful value when dealing with a peaked 
distribution, but none of the actual distributions of |J(x)- 
3(x+r) | has an apparent peak (an example of one of these 
distributions is Figure ITT1 where P(3 ■ 3) is the number 
of haloes in each bin normalized so that the area under 
the curve is unity). A fairer way of dealing with these 
noisy distributions is to fit a straight line and see if there 
is any deviation from randomness. The slope of the best 
fit line indicates if more haloes are aligned parallel or 
orthogonal to each other. 



P(\3(x) ■ 3{x + r)\)= m\3(x) ■ 3{x + r)| 



(15) 



A positive slope (ra) of the best fit line means there are 
more parallel aligned halo pairs, a negative m means they 
are more orthogonal and m = means the haloes have 
random alignment. The values of m that maximized the 
likelihood of fitting a straight line to the distributions 
are shown in the bottom panel of Figure [T2j 

The shape of the plot of the slope (bottom panel of 
Figure fl"2|) is similar to the shape of the plots of the 
conventional statistics. This is expected since they are 
effectively measuring the same thing but in a slightly 
different way. Halo spins are aligned parallel for halo 
separations under 0.3Mpc/h. This alignment has not 
been seen before in simulations because it exists only 
on very small scales which have not before been exam- 
ined. It has howeve r been seen in galax y surveys, for ex- 
ample Galaxy Zoo (|Slosar et al.ll2009() found alignment 
for galaxies closer than 0.5Mpc. The alignment exists 
on the scale of substructure within clusters. If only the 
most massive subhalo in each FOF group is taken (the 
substructure is thrown out), then there is no significant 
alignment at any scale (Figure IT3)) . Here there are no 
haloes at small separations and there is no significant 
alignment at any scale. Only the subhaloes within large 
clusters exhibit any halo-halo spin alignment, although 
it is weak. 
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Fig. 8. — The alignment of dark matter halo velocity with filaments. For all redshifts, haloes are parallel aligned wit h fi laments which 
demonstrates a streaming motion of haloes down bulk flows. Alignment is characterized by the c parameter of Equation 1121 where 9 is the 
angle between halo velocity and filament axis. Lines are colored as in Figure [5] 



4. EVOLUTION OF SPIN PARAMETER 

The spin parameter is a measure of the amount of an- 
gular momentum co ntained in a halo. It was defined in 
IBullock et all (pOOl as. 



A' = 



V2MVR 



(16) 



given the angular momentum J inside a sphere of radius 
R containing mass M, and where V is the halo circular 
velocity at radius R, V 2 = GM/R. 

The distribution of A' over the halos in our sample 
is shown in Figure 1141 It is well fit by a log-normal 
distribution, 



P(A') 



1 



A'V2^cr 



exp 



ln 2 (A7A(,) 
2a 2 



(17) 



The fit was done using a Markov chain Monte Carlo 



maximum likelihood analysis. For all haloes with more 
than 500 particles at z=0 the best fit values are Aq = 



0.02900 
0.02940 



+0.00006 
0.00005' 
+0.00008 
-0.0001 ' 



0.604 



+0.001 
0.002 



and at z=3 An 



a = 0.576 ± 0.002. The distributions 
at both these redshifts over all haloes in the snapshots 
are nearly identical. 

When haloes are binned by mass, the spin parameter 
at high redshift shows a mass dependence while there is 
no mass dependence at z=0, as shown in the left hand 
side of Figure 1151 Here the spin parameter is character- 
ized by the mid point of the log-normal distribution, A . 
The spin parameter over all redshifts is only the same 
for low mass (M < 10 12 ) haloes but there are far more 
low mass than high mass haloes. Since low mass haloes 
dominate, the average distributions over all haloes at the 
different redshifts look the same. At high redshift, there 
is a tendency for the spin parameter to be smaller for 
high mass haloes. 
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Fig. 9. — The alignment of dark matter halo velocity with fila- 
ments on the scale of 2.0Mpc where bulk motions have been sub- 
tracted. Colored lines are for different redshifts as in Figure [5] 
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Fig. 10. — The alignment of dark matter halo velocity with the 
local bulk motion on the scale of 2.0Mpc. Colored lines are for 
different redshifts as in Figure [5] 

This rcdshift dependency can be characterized by a 
power relationship between Aq and mass at each redshift; 

\' <xM at - z \ (18) 

The more negative the value of a, the stronger the cor- 
relation and a = is no correlation at all. The redshift 
dependence of a is shown in Figure [THl The lines for 
haloes with > 500 particles and > 1000 particles overlap 
in Figure [TBI whereas the line for haloes with > 100 par- 
ticles does not. This shows that haloes with more than 
100 particles are susceptible to errors from particles in 
the outer regions and the cut off of only using haloes 
w ith more than 500 part icles is justified. 

iKnebe fc Powerl (|2008[ ) found that mass binning and 
selection criteria for relaxed haloes has almost no effect 
on this correlation. We did find a small effect when a dif- 
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Fig. 11. — An example of a distribution of halo-halo spin align- 
ments. P ^| 3(x) ■ J(x + r)|^ is the number of halo pairs in each bin 

of J ■ J. This example is for haloes that are separated from 0.06 to 
0.1 Mpc/h, which is the second data point from the left in Figure 
1121 The thin line is the actual distribution and the thick line is a 
straight line fit. There is a significant deviation from random spin 
orientations here, shown by the positive slope of the straight line. 
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Fig. 12. — The alignment of neighbouring halo's spins, for haloes 
separate d by distance r. Three differ ent statistics are used: r\ from 
Equation l 131 (top) . rf% from Equation ll4l (middle), and m, the slope 
of the distribution of the halo - halo spin alignment (bottom, an 
example is shown in Figure ITU . The dashed lines are for random 
halo alignments and the shaded regions are the 1<t errors. 

ferent halo catalogue was used. Instead of using all the 
subhaloes, only the most massive subhalo (with more 
than 500 particles) in each friends-of-friends halo was 
used. Most of the mass of the FOF halo is in the most 
massive subhalo so it can be regarded as the background 
halo itself. When substructure is disregarded, we find 
that there is a stronger mass dependency of the spin pa- 
rameter at almost all redshifts (the green line in Figure 
1161 is below the corresponding orange line which includes 
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Fig. 13. — The alignment of neighbouring friends-of-fricnds halo's 
spins. There is no alignment at any scale. 
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Fig. 14. — The distribution of the spin parameter at z=0. The 
histogram is the data and the red (smooth) line is a log-normal fit 
(Equation [TTt where the best fit values are \' = 0.02900±q °0005 > 
„ (i 604+ 001 

all substructure). The spins of subhaloes are greatly af- 
fected by interactions and merger events so may be out 
of equilibrium. 

Mass dependence of the spin parameter a t high redshift 
was first found by IKnebe &: Power! (|2008D . who looked 
at z=l and z=10. When extrapolating the linear trend 
of a(z) with redshift, we predict a much stronger cor- 
relation, a(z — 10) ~ —3 whereas they find a(z = 
10) = —0.059 ± 0.171. Our resu lts agree more closely 
with lMunoz-Cuartas et al.1 (|2011D who found a(z — 2) ~ 
—0.03. For haloes in different environments (blobs, fila- 
ments, sheets and voids), the trends are the same. 

When halo mass is scaled by characteristic mass in 
the right hand side of Figure [15J we find that haloes at 
similar stages of collapse at z=0 and 1 have the same 
spin parameter (the orange and red lines overlap). At 
high redshift, haloes at similar stages of collapse have a 
higher spin parameter (At log M/M* = 3 for example, 
the green (z=3) point lies above the points for z=2 and 
z=l). This may be the result of accretion and merger 
events decreasing the spin of haloes. At z=3, haloes have 
retained much of their initial spin but by z=l, similar 
haloes have experienced accretion that has lowered their 
spin parameter. 

5. SUMMARY AND DISCUSSION 

Using the Millennium N-body simulation, we have 
tracked the evolution of dark matter halo angular mo- 
mentum alignments with the large scale structure, with 
each other and the evolution of the spin parameter. We 



have used the shape of the density field to find filaments 
of 2Mpc in scale in the large scale structure. The align- 
ment between dark matter halo spin and the axis of fila- 
ments was characterized by the shape of the distribution 
of |cos(6>)| where 9 is the angle between the two vec- 
tors. The distribution was fitted to the PDF of Equation 
1121 to find the free parameter c which characterized the 
strength of parallel or orthogonal alignment. 

We found that angular momentum vectors of dark mat- 
ter haloes since z=3 are generally orthogonal to filaments 
but high mass haloes have a stronger orthogonal align- 
ment than low mass haloes. At z=0 the spins of low 
mass haloes have become parallel to filaments, whereas 
high mass haloes keep their orthogonal alignment. 

An interpretation of this is that at early times all halo 
spins were aligned orthogonal to filaments, as TTT pre- 
dicts. High mass haloes especially are well aligned be- 
cause they have had their maximal expansion more re- 
cently and so will have been tidally torqued for longer. 
They usually exist close to clusters where the infall of 
dark matter is almost isotropic and so the nett effect 
from mergers and accretion is minimal. Low mass haloes, 
however, are vulnerable to being disturbed by mergers 
and accretion which is usually assumed to have the ef- 
fect of randomizing the spin orientation. This leaves un- 
explained why low mass haloes at low redshift exhibit a 
parallel alignment with filaments. 

We found that filaments are regions bulk flow. When 
bulk flows are included there is a clear trend for haloes to 
travel parallel to filaments, and high mass haloes travel 
with the best alignment. When bulk flows on the scale 
of the filaments are subtraced, an orthogonal alignment 
to filaments remains, particularly for low mass haloes. 
This shows that entire filaments themselves are moving 
towards attractors and on small scales there is only or- 
thogonal motion. There was also an orthogonal motion 
of low mass haloes with the bulk flow but no alignment 
of high mass haloes out of the bulk flow. 

The motions of haloes relative to the bulk flow could af- 
fect how matter is accreted onto them and the spin orien- 
tation this would cause. Orthogonal motion to the bulk 
flow and filaments by low mass haloes could cause low 
mass haloes to accrete matter preferentially in one direc- 
tion. High mass haloes traveling with the bulk flow would 
experience accretion differently, and this could cause the 
difference in spin orientation. 

Filaments at large smoothing lengths at low redshift 
contain haloes with the best aligned spins and bulk 
motion, while at high redshift it is filaments at small 
smoothing lengths that contain the best aligned haloes. 
This shows that filaments are growing in size over time. 
Because of the nature of the way that the filaments were 
found (using Gaussian smoothing), this enlargement tells 
more about the width of the fila ments rather than the 
length. This is complimentary to iSousbie et al.l (2009) 
where filament length is discussed and it was found that 
there is a general dilation of filaments that began larger 
and a shrinking, fusion and disappearance of the smaller 
filaments. 

We found an alignment only between the spin orienta- 
tion of very close neighbouring haloes. Only at separa- 
tions of less than 0.3Mpc/h do haloes exhibit any mutual 
parallel alignment of their spin axis. The halo finding 
method used in the Millennium simulation has enabled 
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Fig. 15. — The redshift evolution of the spin parameter (Aq). The red line is for z=0, yellow is z=l, blue is z=2 and green is z=3. The 
shaded regions are the la confidence intervals. Left: At high redshift the spin parameter is less and there is a mass dependency. Right: 
Mass bins are normalized by the characteristic mass. 

was done by finding the centre of the log-normal distri- 
bution of the spin parameter. There is a mass dependace 
of the spin parameter at z=3 but not at low redshift and 
the spin parameter is lower overall at high redshift. The 
spin parameter follows a power law with halo mass at 
high redshift but is independent of mass at z=0. 

Future work will bridge the gap between idealistic 
CDM simulations and real galaxy observations. To do 
this we will generate mock galaxy catalogues and use 
only the data that would be available in a real sur- 
vey, see if any alignments of galaxy spin orientations 
could be seen in the universe. This could be used to 
Ian a survey using new multi-object IFU ins truments 
Bland-Hawthorn et aHl2011t iCroom et al.ll2011[ ). 




Fig. 16. — The mass dependence of the spin parameter over red- 
shift. The mas s de pendence, a is the slope of the straight dashed 
lines in Figure [TBI The red (lowest) line includes all haloes with 
more than 100 particles, orange line includes 500 particles and 
blue line includes 1000 particles. The green line is for the halo cat- 
alogue which doesn't include substructure. There is a linear trend 
of stronger mass dependence at higher redshift. 

us to see this small scale alignment. In the Millennium 
simulation, the subfind algorithm was used to identify 
substructure in friends-of-friends groups, and the sub- 
haloes are counted as haloes. This means that align- 
ments between very close haloes can be probed, not just 
alignments between the friends-of-friends groups. 

Lastly, we tracked the evolution of the spin parameter 
from z=3 to now and its dependence on halo mass. This 
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